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ABSTRACT 

We report a supernova discovery in Ks-band images from the NAOS CONICA adaptive optics 
(AO) system on the ESO Very Large Telescope (VLT). The images were obtained as part of a near- 
infrared search for highly-obscured supernovae in the nuclear regions of luminous and ultraluminous 
infrared galaxies. SN 2004ip is located within a circumnuclear starburst at 1.4 arcsec (or 500 pc) 
projected distance from the K-band nucleus of the luminous infrared galaxy IRAS 18293-3413. The 
supernova luminosity and light curve are consistent with a core-collapse event suffering from a host 
galaxy extinction of up to about 40 magnitudes in V-band which is as expected for a circumnuclear 
starburst environment. This is the first supernova to be discovered making use of AO correction and 
demonstrates the potential of the current 8-meter class telescopes equipped with AO in discovering 
supernovae from the innermost nuclear regions of luminous and ultraluminous infrared galaxies. 
Subject headings: supernovae: individual(SN 2004ip) - galaxies: starburst - galaxies:individual(IRAS 
18293-3413) - infrared: galaxies - instrumentation: adaptive optics 



1. INTRODUCTION 

The observed rate at which stars more massive than '^8 
Mq explode as core-collapse supernovae (CCSNe) can 
be used as a direct measure of the current star forma- 
tion rate (SFR), given a reasonable assumption for the 
initial mass function (IMF). In fact, CCSNe are begin- 
ning to be used as probes of the SFR at both low and 
high redshift {e.g. Cappellaro et al. 1999; Dahlen et al. 
2004; Cappellaro et al. 2005) with the aim of providing 
a new independent measurement of the star formation 
history of the Universe. A large fraction of the massive 
star formation at high-z took place in luminous (Ljr > 
10^^ L0) and ultraluminous (Ljr > 10"'^^ Lq) infrared 
(IR) galaxies (henceforth LIRGs and ULIRGs) and their 
high SFRs can be expected to result in CCSN rates a 
couple of orders of magnitude higher than in ordinary 
field galaxies. However, most of the SNe occurring in 
LIRGs and ULIRGs are likely to be obscured by large 
amounts of dust in the nuclear starburst environment 
and have therefore, even in the local Universe, remained 
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undiscovered by optical SN searches. The existence of 
hidden SN factories in the nuclei of LIRGs and ULIRGs 
has already been demonstrated by high-resolution radio 
observations. For example, VLBI observations of the 
nearby ULIRG, Arp 220, have revealed luminous radio 
SNe within its innermost ^150 pc nuclear regions at a 
rate indicating a SFR high enough to power its entire 
IR luminosity (Lonsdale et al. 2006; Parra et al. 2006). 
However, not all CCSNe are likely to become luminous at 
radio wavelengths, and therefore radio searches for SNe 
can only provide a lower limit for the real SN rates. 

In the near-IR K-band the extinction is strongly re- 
duced while the sensitivity of ground-based observations 
is still good, making searches for SNe in starburst galax- 
ies, LIRGs and ULIRGs more promising (Van Buren et 
al. 1994a; Grossan et al. 1999; Mattila & Meikle 2001 
henceforth MMOl; Maiohno et al. 2002; Mannucci et al. 
2003; Mattila et al. 2004). Furthermore, the recent in- 
troduction of adaptive optics (AO) for 8-meter class tele- 
scopes now enables near-IR searches for SNe at a spatial 
resolution comparable to some of the radio studies. In 
this letter we report a SN discovery in the circumnuclear 
starburst of IRAS 18293-3413 which is a LIRG at a red- 
shift of z = 0.0182 (Strauss et al. 1992) corresponding to 
a distance of 79 Mpc (Hq = 70 km s-iMpc-\ Oa=0.7 
and riM=0.3). The SN was detected as a result of our 
Ks-band search of highly obscured CCSNe in a sample of 
15 nearby LIRGs and ULIRGs using the NAOS CONICA 
(NACO) AO system on the ESO Very Large Telescope 
(VLT). 

2. OBSERVATIONS AND RESULTS 

2.1. NACO observations and data reduction 

IRAS 18293-3413 was observed for the first time with 
VLT/NACO (Rousset et al. 2003; Lenzen et al. 2003) in 
service mode on 2004 May 4.3 UT in the Ks-band with 
the S27 camera (0.027 arcsec/pixel) using the autojitter 
imaging sequence (see Table 1 for details on the observ- 
ing parameters and conditions). The AO correction was 
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TABLE 1 

Log of VLT/NACO observations of IRAS 18293-3413. 



Date (UT) 
2004 



Band 



Integration 



Seeing Airmass SR 

(arcsec) (%) 

May 4.3 24 X 90 0.47-0.79 1.08-1.19 20-30 (22) 

Sept. 13.1 Ks 23 x 90 0.65-0.90 1.15-1.28 15-20 (21) 

Sept. 15.0 Ks 63 X 90 0.70-1.25 1.01-1.15 5-20 (41) 

Sept. 27.1 H 68 x 90 0.53-1.15 1.11-1.70 5-15 (38) 

Sept. 27.1 Ks 15 x 90 0.69-0.81 1.81-2.13 5-10 (6) 

Note. — The strehl ratios (SR) of frames used for the combined 
images are listed in col. 6 with the number of frames in brackets. 

performed using the visual wavcfront sensor (VIS-WFS) 
with a natural guide star of my "^15, 25" northeast from 
the science target. The observations were repeated on 
2004 Sept. 13.1 UT to allow SN detection by subtract- 
ing the images from the two epochs. The data obtained 
on Sept. 13 were immediately reduced and searched for 
any obvious SNe (for details see below and Section 2.2). 
This analysis produced a clear detection of a new point 
source in the Sept. 13 data. We therefore triggered the 
Target of Opportunity (ToO) part of the programme to 
obtain a deep repeat observation in the Ks-band on Sept. 
15.0 UT, two days after the discovery image. About 
two weeks later, on Sept. 27.1 UT we triggered VLT 
ToO again, this time to obtain both light curve and H-K 
colour information for the SN. The SN was reported on 
the 'CBAT Possible Supernovae Page'-'^ immediately af- 
ter the discovery and was later assigned the designation 
SN 2004ip (Mattila et al. 2007). 

The NACO data were reduced using the IRAF DIM- 
SUM package and IDL. The jittered on-source frames 
were median combined to form a sky frame. The sky 
subtracted images were de-dithered making use of the 
centroid coordinates of a bright field star visible in all 
the frames, and the de-dithered frames were median com- 
bined excluding frames with the lowest strehl ratios (SR) 
(see Table 1). 

2.2. Supernova detection and photometry 

The images from Sept. 13, Sept. 15, and Sept. 27 were 
rotated and shifted to align them to the reference im- 
age (from May 4) using bright isolated stars around the 
galaxy. We used a slightly modified version of the Opti- 
mal Image Subtraction method (Alard & Lupton 1998, 
Alard 2000) as implemented in ISIS 2.2 to derive a con- 
volution kernel to match the better seeing image to the 
image with the poorer seeing and to match the intensi- 
ties and background differences between the frames. In 
MMOl we demonstrated the feasibility of this method 
in nuclear SN detection. The subtracted images show a 
clear point source (with its peak at the ~20cr and ~7cr 
level in Sept. 13-15 and Sept. 27 epochs, respectively) 
at 1.14" east and 0.78" north from the Ks-band nucleus 
of the galaxy (see Fig. 1). However, without a suitable 
reference image for image subtraction we were unable 
to detect the SN in our H-band data. Using the world 
coordinate system from a 2MASS image of IRAS 18293- 
3413 we obtained R.A. = 18h32m41s.26 and Decl. = 
-34''11'26.7" (equinox 2000.0) for the SN with an esti- 
mated precision of ±0.4". 

We adopted the photometric zeropoints provided by 
ESO for our nights of observations and estimate them 

1 http://cfa-www.harvard.edu/iau/CBAT_PSN.html 



to be accurate to ±0.03 magnitudes in comparison to 
2MASS data. Since the SR changes as a function of the 
distance from the AO guide star, to measure the SN pho- 
tometry with the associated uncertainties, we used two 

-isolated bright stars lying, respectively, 12" and 38" from 
the guide star (7.1"N, 18.3"E and 5.8"S, 11.5"W of the 
galaxy K-band nucleus) to provide a point spread func- 
tion (PSF) with a higher and lower SR compared to the 

-actual SN. The two PSFs were scaled to correspond to 
magnitudes ranging between 18.5 and 19.1 in steps of 0.1. 
Sources with each magnitude and PSF were then simu- 
lated in the images at twenty different locations evenly 
distributed around the galaxy nucleus at the same offset 
as the real SN. The sources were added to the images one 
at a time before the images were aligned, matched and 
subtracted. The magnitudes for both the real and sim- 
ulated SNe were meastired from each subtracted image 
using aperture photometry with a 0.15" aperture radius, 
and a sky annulus between 0.20" and 0.30". We then 
estimated a magnitude for the SN by comparing its pho- 
tometry with the mean of the counts measured for the 
20 simulated sources with each magnitude and PSF. This 
yielded SN magnitudes of m(Ks) = 18.74 and 18.64 for 
Sept. 13 using the two different PSFs, respectively. Sim- 
ilarly, we obtained m(Ks) = 18.66 and 18.72 for Sept. 
15. We adopted the standard deviation (a) of the 20 
measurements as the uncertainty in our SN photometry, 
yielding a = 0.12 and 0.09 for the two epochs, respec- 
tively. Measuring with a 50% larger aperture and sky 
annulus gave very similar magnitudes as expected since 
the photometry is performed relative to the simulated 
sources, only the uncertainties were increased. During 
the Sept. 27 Kg-band observation the airmass was high, 
resulting in a low SR and therefore wider PSF for the im- 
ages. For this epoch stronger residuals were also present 
in the subtracted images and we therefore excluded ten 
of the twenty simulated sources which were affected by 
much stronger residuals than the SN. For the Sept. 27 
data we used an aperture radius of 0.30" , and a sky an- 
nulus between 0.40" and 0.50". This yielded m(Ks) = 
18.82 and m(Ks) = 18.99 with a = 0.14 using the two 
PSFs, respectively. We adopt the average of the magni- 
tudes obtained with the two different PSFs as the best 
estimate for the SN magnitude. The uncertainties were 
obtained as a quadrature sum of the a value and the 
error in the zeropoint. 

2.3. Supernova type and extinction 

Without a spectroscopic observation of the SN its con- 
firmation and possible typing have to be based on pho- 
tometry and other available information. An asteroid 
origin for the source was immediately ruled out since the 
images separated by 2 days show the SN in exactly the 
same position (within ±0.02"). Furthermore, the SN is 
still at precisely the same location in the 3rd epoch im- 
ages, obtained 2 weeks from the discovery. Below we 
model the near- to far-IR SED of IRAS 18293-3413 and 
find a negligible contribution from an active galactic nu- 
cleus (AGN) to its Imninosity. Also, the presence of an 
AGN is not favoured by the X-ray results of Risaliti et al. 
(2000) and the fact that IRAS 18293-3413 shows an 'H 
II galaxy' spectrum in the optical (Veilleux et al. 1995). 
Therefore, an AGN flare can be ruled out as origin for 
the source. A variable foreground star also provides a 
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Fig. 1. A 16.5" x 11.5" subsection of the NACO Ks-band image of IRAS 18293-3413 obtained on 2004 Sept. 15.0 UT with North up 
and East to tlic left (A). The result of image matching and subtraction of the NACO Kg-band reference image (May 4.3) from the images 
obtained on Sept. 13.1 (B) and Sept. 15.0 (C). SN 2004ip (indicated by tick-marks) is clearly detected in the subtracted images at 1.14" 
east and 0.78" north of (or 500 pc projected distance) from the K-band nucleus. Simulated sources of m(Kg) = 18.7 are shown southeast 
and northwest of the SN in the two subtracted images, respectively. A strong negative residual is visible at the location of the galaxy 
nucleus in the subtracted images. All the images are shown with an inverted brightness scale. 



very improbable alternative explanation given the small 
field of view (42"x42", i.e., 1/(3x10**) of the whole sky) 
of our combined NACO frames. For example, the likeli- 
hood of catching a foreground nova would be extremely 
small given the Galactic nova rate of ~35 yr~^ (Shafter 
1997). Furthermore, a classical nova in IRAS 18293-3413 
would be several magnitudes fainter than the detected 
SN (Hachisu & Kato 2006). We, therefore, conclude that 
anything other than a SN origin for our source appears 
extremely unlikely. 

To estimate the SFR for IRAS 18293-3413 we compiled 
its 1-100 micron archival (2MASS and IRAS) photometry 
and fitted the resulting SEDs (sec Fig. 2) both with the 
starburst models of Efstathiou et al. (2000), and the cir- 
rus model of Efstathiou & Rowan- Robinson (2003). The 
combination of the starburst and cirrus model provides a 
very good fit to the near- to far-IR SED and is consistent 
with a negligible contribution from an AGN. The model 
assumes an exponentially decaying starburst with a time 
constant of 40 Myrs and an age of 70 Myrs, and an un- 
derlying 12.5 Gyr old stellar population. The cirrus Ay 
is assumed to be 1.7 and the initial Ay of the molecular 
clouds that constitute the starburst is 100. The SFR at 
the peak is 285 Mq yr^^, and 135 Mq yr~^ when aver- 
aged over the duration of the starburst. We estimate a 
rest frame 8-1000 /im luminosity of 5.8 x 10^^ Lq and a 
bolometric luminosity of 6.5 x 10^^ Lq for IRAS 18293- 
3413. 

Adopting the average SFR and assuming CCSN pro- 
genitors between 8 and 50 and a Salpeter IMF with 
cut-offs at 0.1 and 125M0 (see also MMOl), we obtain 
a CCSN rate of -1 yr"! for IRAS 18293-3413. We note 
that this CCSN rate also agrees with the model predic- 
tions of Genzel et al. (1998) for an exponentially decay- 
ing starburst with similar time constant, age and bolo- 
metric luminosity. Recent observational evidence shows 
that the Type la SN (SN la) rate is likely to depend on 
both the stellar mass and the mean SFR of their host 
galaxies (sec Sullivan et al. 2006 and references therein) . 
The rate of SNe la associated with the old stellar popula- 
tion is not expected to be enhanced in IRAS 18293-3413 
compared with normal field galaxies but the rate asso- 
ciated with the recent star formation could be higher. 
This rate has becni found to depend on the host galaxy 
SFR averaged over the last 500 Myrs (SFRave) accord- 
ing to ~4 X 10"^ SNe yr-i / SFRave (Sullivan et al. 
2006). Using 135 M© yr~-^ as an upper limit for SFRave 
we obtain an upper limit for its SN la rate of ~0.05 yr~^ 
which is ~20x smaller than the estimated CCSN rate. 
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Fig. 2.— The IR SED of IRAS 18293-3413 fitted with the 
starburst models of Efstathiou et al. (2000) and the cirrus model 
of Efstathiou and Rowan- Robinson (2003). No AGN contribution 
was needed to explain the SED. 

We therefore conclude that SN 2004ip was very likely a 
core-collapse event. 

After correcting for a Galactic extinction of Ay — 0.47 
(Schlegel et al. 1998) the absolute magnitude of the SN 
at the time of the discovery becomes Mks = —15.85 (be- 
fore correcting for the host galaxy extinction). In Fig. 3, 
we compare the photometry of the SN with template K- 
band light curves for ordinary and slowly-declining CC- 
SNe from MMOl. These templates comprise 2 and 3 
linear components, respectively, and are based on the 
light curves of 11 (type II and Ib/c) and two (the type 
IIL SN 1979C and the type Iln SN 1998S) CCSNe. The 
slow decliners are significantly more luminous than the 
ordinary events, even at early times, therefore requir- 
ing a larger host galaxy extinction to match the obser- 
vations. Assuming that the SN was discovered at the 
K-band maximum light, we need to dim the ordinary 
and slowly declining templates by Ak = 2.75 and 4.15 
magnitudes, respectively, to match the observations. To 
estimate the lowest plausible extinction towards the SN 
we sought the latest possible SN epoch of discovery to 
match the templates. The time interval between the ref- 
erence and discovery images was ^130 days and the tem- 
plate light curves of MMOl have rise times of ~25 days 
to the K-band maximum light. The latest SN epoch of 
discovery assuming that the SN exploded soon after the 
observation of the reference image is therefore about 100 
days past the K-band maximum light. We found that 
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Fig. 3. — Absolute magnitudes of SN 2004ip compared with 
template K-band light curves (MMOl) for ordinary (solid line) and 
slowly-declining (dotted line) CCSNe. The SN epoch of discov- 
ery was assumed to be at (K-band) maximum light (A), 100 days 
past the maximum light (B). The host galaxy extinctions used for 
dimming the light curves are shown in the figures. 

the ordinary and slowly-declining templates can match 
the observations at this epoch if Ak = 0.55 and 2.50 
magnitudes, respectively. Adopting the extinction law of 
Rieke & Lebofsky (1985) the range of likely host galaxy 
extinction towards SN 2004ip is therefore between about 
5 and 40 magnitudes in Ay. 

3. DISCUSSION 

The number of CCSNe discovered in the infrared (Van 
Buren et al. 1994b; Maiolino et al. 2002; Mattila et 
al. 2005a, b) is still rather small. However, these few 
events can have an important impact for the SN statistics 
when estimating the complete local CCSN rate includ- 
ing also the SNe in the optically obscured parts of the 
galaxies. Detecting CCSNe directly in starburst galaxies 
over a range of IR luminosities (Lir = 10^°-10^^ Lq) is 
also vital for interpreting the results of the high-z CCSN 
searches since a large fraction of the massive star forma- 
tion at high-z took place in IR luminous galaxies. 

SN 2004ip is the first reported supernova discovered 
using adaptive optics. This discovery demonstrates the 
great potential of the current 8-meter class telescopes 
equipped with AO correction in discovering SNe from 
the innermost nuclear regions of LIRGs and ULIRGs. In 
addition, the next generation wide-field AO systems can 



be used for expanding near-IR SN searches to the high-z 
Universe. Near-IR-discovered SNe in nuclear starbursts 
are not only important for SN rate estimates, they can 
also be used to estimate the amount and distribution of 
extinction within their host galaxies, and follow-up ob- 
servations can provide us with a better understanding 
of the behaviour of SNe within the dusty, high-density 
starburst environment. The extinction of up to 40 mag- 
nitudes in Ay estimated for SN 2004ip is consistent with 
those derived towards the radio SNRs in the nuclear (cen- 
tral 500 pc) regions of M 82 (MMOl). We note that 
Infrared Space Observatory (ISO) studies (Genzel et al. 
1998) indicate much higher extinctions of Ay (mixed) ~ 
50-1000 towards the highly-excited mid-IR line-emitting 
innermost nuclear regions in LIRGs and ULIRGs. How- 
ever, SN 2004ip is located in a circumnuclcar starburst 
at a 500 pc distance from the Kg-band nucleus where the 
conditions could be expected to be similar to the M 82 
nuclear environment. 

Many of the SNe exploding in the nuclear starburst 
environment can be expected to produce very luminous 
radio SNe as a result of the interaction of the SN ejecta 
with a dense surrounding medium (Chevalier and Frans- 
son 2001). In fact, the radio discovered SN 2000ft (Col- 
ina et al. 2001; Alberdi et al. 2006) in a circumnu- 
clcar starburst ring (600 pc from the nucleus) of its LIRG 
host NGC 7469, showed a luminosity comparable to the 
brightest known radio SNe. Radio follow-up observation 
of SN 2004ip and other SNe discovered in LIRGs and 
ULIRGs in the future can therefore provide important 
information on their progenitor stars and on the condi- 
tions in the surrounding circumstellar medium. 
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0406). We thank the ESO staff for carrying-out the ser- 
vice observations and C. Alard and R. Joseph for useful 
discussions. 
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